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I.  INTRODUCTION 


This  report  presents  an  analysis  of  viscously  damped  tracking  aids. 
The  problem  begins  when  a human  needs  to  manually  point  something  like 
a weapon  or  a camera.  The  pointing  accuracy  can  be  improved  if  the 
object  is  supported  by  something  like  a tripod  and  if  viscous  damping 
is  added  to  smooth  out  the  tracking.  This  report  applies  classical  con- 
trol theory  to  the  problem.  It  supports  the  argument  with  a series  of 
simulations  with  a human  in  the  loop.  It  gives  some  preliminary  design 
guidance  for  stiffness  and  damping.  Finally,  it  discusses  some  general 
applications  and  requirements  for  future  research. 

Viscously  damped  tracking  aids  have  been  used  since  1962.  At  that 
time  the  US  Army  Human  Engineering  Laboratories  started  using  them  for 
directing  anti-tank  guided  missiles.1’2’3  Currently  there  are  approxi- 
mately four  Army  applications  and  at  least  one  commercial  camera  tripod 
that  use  viscous  damping.  They  have  been  designed  with  the  general  notion 
that  there  should  be  a lot  of  damping  so  as  to  mask  the  effect  of  inertia; 
that  the  legs  should  be  as  stiff  as  possible  to  prevent  windup  and  oscil- 
lations; and  that  the  coulomb  friction  should  be  as  low  as  possible. 

The  amount  of  inertia  has  always  been  a fixed  parameter  determined  by  the 
device  being  pointed.  Several  iterations  of  build  and  test  have  achieved 
very  fine  performance  indeed.  For  example,  a nominal  set  of  specifica- 
tions might  be  as  follows: 


Inertia 

J = 

1 slug-ft^ 

Damping 

B = 

100  ft- lb/rad/sec 

Stiffness 

K = 

10000  ft- lb/rad 

Accuracy 

RMS  = 

100  microrad  (low  tracking  rates) 

Until  now  no  one  has  performed  the  classical-control  analysis  to  make 
determinations  of  design  parameters. 


II.  SYSTEM  MODEL 

The  first  step  was  to  develop  a model  of  the  control  loop.  The 
lumped  parameter  model,  block  diagram,  and  transfer  functions  shown  in 


*R.  T.  Gschwind,  "Gunner  Tracking  Behavior  as  a Function  of  Three 
Different  Control  Systems,"  Human  Engineering  Laboratories  Technical 
Memorandum  2-63,  January  1963.  */0V  OSS’ 

2 

R.  T.  Gschwind,  "Gunners'  Aiming  Errors  in  Antitank  Weapons,"  Human 
Engineering  Laboratories  Technical  Memorandum  5-64,  Feb  64.  /13>-  3iT/  3*7 O. 

3 

F.  N.  Newcomb,  "Viscous-Damping  Mechanisms  as  Applied  to  4-Inch  Rocket 
Launcher  Mount,"  Human  Engineering  Laboratories  Technical  Note,  Feb  63. 
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f PRECEDING  PASS  BLANK- NOT  FILMED 


Figure  1 are  the  result  of  the  following  rationale  and  simplifying 
assumptions: 

a.  It  is  a one-dimensional  tracker. 

b.  The  human  applies  torque  directly  to  the  load  through  a fixed 
handle;  therefore,  there  is  no  gear  reduction  or  backlash. 

c.  All  inertia,  J,  is  located  in  the  rotating  mass. 

d.  All  compliance,  lCp,  is  located  between  the  damper  and  the 
ground. 

e.  The  viscous  damping,  B,  is  linear  with  angular  rate. 

f.  There  is  no  coulomb  friction. 

g.  An  optical  sight  provides  a visual  error  signal  to  the  human 
representing  the  angular  error  between  the  device  angle,  and  the  ref- 
erence (target)  location,  r. 

h.  The  sight  provides  enough  magnification  that  the  eyeball  thresh- 
old is  not  reached. 

i.  Whatever  the  sight  magnification  is,  it  is  included  in  the 
human  gain  K^. 

j . The  human  output  torque  is  linear  with  respect  to  error  regard- 
less of  magnitude. 

A human  transfer  function  was  chosen  from  "Man-Machine  System"  by 
Sheridan  and  Ferrell.  Actually,  Sheridan  selected  it.  He  believed 
for  a simple  task  such  as  a heavily  damped  tripod  with  a low  frequency 
reference  task  that  a pure  lag,  t = 0.3  seconds,  and  a gain,  1^,  would 

be  a good  approximation  of  a human;  furthermore,  the  Pad§  Approximation 
can  be  used  for  such  a low  frequency  problem. 

Yh  = 1^  exp  (-ts)  ■ (1  - J S) /Cl  + \ S) 

The  transfer  function  resulting  from  these  assumptions  is  shown  in 
Figure  1.  It  was  used  in  the  several  analyses  that  follow. 


III.  ROUTH  STABILITY  CRITERION 

Routh  Stability  Criterion  was  applied  to  the  characteristic  equa- 
tion to  find  the  upper  bound  on  human  gain,  K^,  as  a function  of  the 

tripod  parameters  J,  lCj,,  and  B.  This  analysis  presumes  that  the  human 
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can  increase  his  gain  almost  without  limit  until  he  is  satisfied  with 
the  system  performance  or  the  system  goes  unstable.  Routh's  method  was 
an  easy  way  to  find  the  stability  limit  for  a wide  range  of  parametric 
values.  Actually,  some  scaling  was  found  to  help  in  simplifying  the 
presentation  of  these  results. 

a.  The  velocity  error  coefficient,  K^/B, was  used  rather  than  1^. 

b.  The  closed  loop  zero.-lCjVB,  was  used  because  several  analyses 
show  it  to  have  a sensitive  and  consistent  effect  on  performance. 

c.  The  final  term  was  chosen  as  /lCp/J  because  it  gets  the  job 

done;  i.e.,  it  includes  the  effect  of  J,  it  reduces  the  range  of  K^,,  it 

appears  in  the  transfer  function,  and  it  approximates  the  imaginary 
dimension  of  a pair  of  closed  loop  poles. 

Figure  2 shows  the  maximum  allowable  values  of  K^/B  as  determined 
from  Routh's  Criterion.  First,  it  shows  that  lC^/B  needs  to  be  greater 
than  6.6,  which  incidentally  corresponds  to  the  open  loop  pole-2/x. 

Second,  it  shows  that  maximum  K^/B  is  achieved  with  lCp/B  equal  to  14  and 

as  large  as  possible.  Of  course,  maximum  K^j/B  is  not  neces- 
sarily the  same  as  minimum  tracking  error  as  will  be  seen  in  the  next 
two  analyses. 


with  ZlCp/J 


IV.  INTEGRAL  SQUARE-ERROR  (ISE)  CRITERION 

The  ISE  Criterion  was  used  to  see  how  the  choice  of  system  parameter 
values  can  affect  tracking  accuracy.  The  ISE  was  selected  because  it 
can  be  calculated  from  the  transfer  function,  thereby  permitting  easy 
variation  of  parameter  values.  The  solution  of  the  integral  can  be 
found  in  Analytical  Design  of  Linear  Feedback  Systems,  by  Newton,  1957. 
For  convenience  the  computational  algorithm  is  presented  on  a separate 
page  (Table  I). 

A FORTRAN  program  was  written  which  varied  human  gain  until  a mini- 
mum value  of  ISE  was  found  for  each  combination  of  tripod  parameter 
values.  These  results  are  plotted  in  Figure  3.  The  ISE  curves  show  a 
ratio  of  ICp/B  equal  to  10  to  be  best  for  high  values  v^ICp/j.  They  show 

diminishing  benefits  for  increases  in  /k^7J  much  beyond  25  (note  the 

ordinate  does  not  go  to  zero) . The  values  for  K^/B  corresponding  to 

minimum  ISE  are  not  shown;  however,  a value  of  3.0  was  near  optimum  for 
nearly  all  of  the  better  combinations. 
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9 


gure  3.  I SI: 


where  e(t)  = error  as  a function  of  time 

C(s)/D(s)  * frequency  domain  representation  of  e(t). 

2 3 

C(s)  = cq  + CjS  + c2s  + c^s 

D(s)  = dQ  + dj^s  + d2s2  + d3s3  ♦ d4s4 

Co  = 1(T 

cx  = (J/B  + t/2)  ICj. 
c2  = (1  + t/2  • Kj./BJJ 
Cj  * J • t/2 

dQ  = KhVb 

dl  = Kt  ♦ Kn  - Kh(Kt/B)t/2) 
d2  = CJ  Kf/ B ♦ ICj.t/2  - Ktf/2) 

d3  = c2 
d4  = c3 

ISE  = {c2  (-d2d3  ♦ dod1d2)  + (c2  - 2Clc3)  dod1d4 

* (C?  ' Vi>  dod3d4  * Co'-dld4  * d4d3d2>} 
{2dod4  <-dod3  - dld4  * W} 
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The  ISE  is  not  a very  selective  measure  of  performance  as  compared 
to  the  ITAE  considered  next.  These  results  do,  however,  give  a good 
general  picture  of  how  performance  is  affected  by  the  parameter  varia- 
tions, and  they  give  a good  starting  point  for  the  ITAE  analysis,  namely 

lCj./  B = 10 

= 30 

1^/B  = 3.0 


V.  INTEGRAL-OF-TIME-MULTIPLIED  ABSOLUTE  ERROR  (ITAE)  CRITERION 

The  ITAE  Criterion  is  more  selective  than  the  ISE.  It  will  select 
a system  that  damps  out  faster  than  the  ISE  because  it  multiplies  the 
error  by  time.  The  bad  feature  of  the  ITAE  is  that  it  must  be  deter- 
mined by  measurement  of  the  system  response  unlike  the  ISE  and  Routh's 
Criterion  which  may  be  calculated  from  the  transfer  function.  The  ITAE 
was  determined  from  the  closed  loop  response  to  a unit  step  as  calcula- 
ted by  DYSYS,  a Fourth  Order  Runge-Kutta  Routine.  A calculation  was 
made  every  tenth  second  for  four  seconds  after  the  step.  The  absolute 
error  was  multiplied  by  time  and  by  the  time  increment  and  summed  up  to 
get  the  ITAE.  Figure  4 shows  two  step-responses  which  indicate  that  the 
four-second  duration  approached  infinity  for  purposes  of  this  problem. 

The  scaled  parameters  were  varied  one  at  a time  until  a minimum 
ITAE  was  achieved.  This  minimum  value  was  used  while  another  parameter 
was  varied  until  a new  minimum  was  found,  and  so  on.  The  data  collected 
during  this  procedure  are  listed  in  Table  II.  Many  of  these  data  are 
plotted  in  Figure  5.  The  bottom  curve  shows  the  ITAE  to  be  very  selec- 
tive in  that  there  is  a clearly  defined  minimum  at  1Cj./B  equal  to  eight. 

Further  changes  in  the  other  two  parameters  could  produce  no  more  reduc- 
tion. There  is  always  the  possibility  that  the  coarseness  of  this  anal- 
ysis could  overlook  the  true  minimum,  but  for  now  the  ITAE  Criterion 
indicates  an  optimum  design: 

JCj,/B  = 8.0 
vic^TT  = 30 

Remember  that  this  optimum  is  for  a step  input  and  a particular  human 
transfer  function.  Additional  support  and  verification  are  required 
before  concluding  it  is  the  optimum  design. 


VI.  ROOT  LOCUS  METHOD 

Root  locus  plots  help  to  understand  why  the  tripod  legs  should  have 
some  compliance  and  why  there  is  an  optimum  value  of  ICp/B.  Figure  6 was 
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se  to  a Step  Input 


TABLE  II.  List  of  Conditions  for  ITAE 


VB 

ICp/B 

/k^Tj 

ITAE 

3.0 

10 

15 

.5093 

3.0 

10 

22 

.3247 

3.0 

10 

30 

.3140 

3.0 

10 

40 

00 

2.5 

10 

22 

.4006 

3.5 

10 

22 

.3306 

2J 5 

10 

30 

.3910 

3.5 

10 

30 

.2782 

3.0 

8 

22 

.3207 

3.0 

12 

22 

.4574 

3.0 

8 

30 

.2327 

3.0 

12 

30 

.4037 

3.5 

8 

22 

.5826 

3.5 

8 

30 

.1781 

4.0 

8 

30 

.2243 

3.5 

7 

30 

.2611 

3.5 

9 

30 

.2222 

3.5 

8 

35 

.1780 
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drawn  to  help  motivate  this  understanding  even  though  the  graphs  are  not 
to  exact  scale. 

Figure  6(a)  shows  the  root  locus  for  a system  with  infinite  stiff- 
ness. The  pole  at  the  origin  is  a free  integration  unaffected  by  any- 
thing we  do.  The  pole  and  zero  at  -6.6  and  +6.6  are  the  Pad6  approxima- 
tion of  a man  assumed  to  be  independent  of  tripod  parameters.  The  pole 
at  -25  is  determined  by  the  damping  B.  This  plot  shows  that  the  locus 
near  the  origin  is  not  going  to  be  affected  much  by  the  value  of  B as 
long  as  it  is  greater  than  ten.  The  human  in  the  loop  will  tend  to 
adjust  his  gain  until  the  closed  loop  root  nearest  the  origin  is  near 
the  dashed  line  indicating  a damping  ratio  of  0.7.  In  general,  the  sys- 
tem response  is  improved  if  this  root  can  be  moved  further  away  from 
the  origin  but  still  near  the  dashed  line  and  of  course  high  gain  is 
desirable. 


Figure  6(b)  shows  the  locus  of  roots  for  a system  with  finite  stiff- 
ness. An  open  loop  zero  is  located  at  -K^/B  and  a pair  of  open  loop 


poles 


are  located  at  -1Cj./2B 


For  the  condition  shown 


here,  K_  has  been  reduced  just  enough  to  make  the  pair  of  open  loop  poles 

1 2 
at  -ICp/2B  become  imaginary;  i.e.,  ICp  is  slightly  less  than  (2B)  /J.  If 

> (2B)2/J,  this  pair  of  open  loop  poles  would  be  real.  The  closed 


loop  root  near  the  origin  will  still  be  near  the  intersection  with  the 
dashed  line  and,  therefore,  it  will  not  be  affected  significantly  by 
small  changes  in  ICp  or  B. 


Figure  6(c)  shows  what  happens  after  K^.  has  been  reduced  even  more. 

The  segment  of  locus  crossing  the  imaginary  axis  (the  stability  limit) 
is  now  controlled  by  the  open  loop  poles  located  at  -1CJ./2B  + j /ICp/J 

(approximate  for  reasonably  large  values  of  say  greater  than  15) . 

The  closed  loop  root  nearest  the  origin  and  along  the  dashed  line  is 
influenced  by  the  location  of  the  zero  at  -ICj./B.  Now  we  can  really 

affect  the  system  performance  by  careful  choice  of  these  parameters. 

The  plot  in  Figure  6(c)  is  getting  near  optimum  because  the  intersection 
of  the  locus  with  the  dashed  line  is  able  to  move  further  away  from  the 
origin.  In  fact  since  a human  can  alter  his  gain  somewhat,  the  actual 
operating  point  might  be  anywhere  in  the  region  of  near  tangency. 


Figure  7 shows  the  root  locus  for  the  so-called  ITAE  optimum 
design.  By  increasing  B the  zero  at  -iCp/B  was  moved  closer  to  the 

origin  and  the  locus  "circle"  became  smaller,  causing  the  system  to  be 
slightly  more  damped  than  the  system  in  Figure  6(c).  The  closed  loop 
root  determined  by  the  ITAE  had  K„/B  = 3.5.  This  point  would  fall  about 
halfway  between  the  closed  loop  roots  shown.  This  system  has  good 
damping  and  good  frequency  response. 
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VII.  SIMULATION 


The  simulation  task  attempted  to  validate  the  transfer  function 
assumed  for  a human  operator  and  to  verify  the  trends  found  in  the  anal- 
yses. The  effort  consisted  of  developing  an  analog  computer  simulation 
of  the  tripod  dynamics,  writing  a hybrid  program  to  operate  the  analog 
and  analyze  the  results,  and  interfacing  a strain  gage  type  force  con- 
troller with  the  analog  computer.  Figure  8 shows  the  physical  layout 
and  the  essential  part  of  the  analog  circuit. 

The  experimental  scheme  required  a subject  to  observe  an  error  sig- 
nal on  a CRT.  He  attempted  to  reduce  the  error,  i.e.,  return  a dot  to 
the  center  of  the  display,  by  applying  a force  to  the  strain  gage  con- 
trol stick.  The  gage  output  signal  entered  the  analog  computer  as  a 
torque  to  drive  the  tripod.  The  tripod  output  angle,  6,  was  subtracted 
from  the  reference  (filtered  white  noise  source)  and  presented  on  the 
CRT. 


The  reference  task  was  selected  to  meet  two  criteria.  First,  the 
task  must  not  be  predictable  by  the  subject  as  a simple  sine  wave  would 
be.  Second,  it  must  distribute  the  tracking  energy  evenly  over  the  fre- 
quencies of  interest  so  that  the  transfer  function  could  be  calculated. 
These  criteria  were  met  by  using  a gaussian  random  variable  for  the 
target  acceleration.  A new  value  was  introduced  every  time  step  (mean 

0. 0,  standard  deviation  0.1).  This  value  was  passed  through  a second 
order  filter  to  keep  the  energy  requirement  for  input  torque  constant 
across  the  frequencies  of  interest.  It  was  also  necessary  to  reduce 
the  low  frequency  drift  to  prevent  overloading  the  filter  integrator 
as  shown  in  Figure  8.  The  filter  transfer  function  then  became 
1/S(S  +1). 

The  data  analysis  was  developed  to  calculate  the  human  gain  and 
phase  lag  every  radian  per  second  from  zero  to  ten.  The  desired 
resolution  of  one  radian  per  second  determined  that  there  be  approxi- 
mately one  minute  of  continuous  data  to  sample  ten  cycles.  The  maximum 
frequency  of  ten  radians  per  second  required  that  there  be  at  least 
twenty  samples  per  second  to  prevent  aliasing,  i.e.,  about  ten  times 
the  frequency.  The  Fast  Fourier  Transform  subroutine  (FFT)  required 

2N  samples.  With  N = 9 there  are  1024  samples  equivalent  to  51  seconds 
of  data  at  20  samples  per  second.  This  sample  size  gives  90  percent 
confidence  of  accuracy  within  25  percent  when  there  are  ten  frequency 
bands  of  interest.  Therefore,  each  trial  on  the  computer  consisted  of 
a 51-second  data  run  plus  a two-second  warmup  period  to  achieve  steady 
state.  Every  time  step  (0.05  seconds)  the  error  and  the  torque  were 
A/D  converted  and  stored  and  a new  Gaussian  sample  was  drawn  and  delivered, 

1. e.,  D/A  converted. 

The  FFT  subroutine  was  called  to  process  both  the  error  and  torque 
data.  In  general  the  ratio  of  the  output  (torque)  to  input  (error)  will 
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yield  the  gain  and  phase  as  a function  of  frequency  from  which  the  trans- 
fer function  can  be  constructed.  In  practice  it  was  necessary  to  combine 
the  vector  pairs  from  eight  frequency  bands  to  obtain  a one  radian  per 
second  band.  It  was  also  necessary  to  rotate  the  vector  pairs  to  a 
common  reference  (zero)  to  avoid  going  past  360°  when  calculating  phase. 
The  rotation  was  performed  first,  then  the  vector  addition,  then  the 
division.  The  remaining  vector  contains  the  gain  and  phase.  More 
averaging  was  achieved  by  running  two  repetitions  for  each  of  two  subjects 
on  each  test  condition.  The  averaged  results  for  the  four  tests  of 
each  condition  are  presented  in  Figures  9 and  10. 

* The  test  conditions  were  chosen  to  give  a wide  variation  in  param- 

eters based  on  the  prior  analysis  and  also  evaluate  some  specific  near 
optimum  designs.  The  four-test  condition  had  parameter  values  as  shown 
below: 


Condition 

J 

B 

*r 

JCjVB 

/iyj 

I 

1 

50 

200 

4 

14 

II 

1 

50 

500 

10 

22 

III 

1 

50 

00 

oo 

OO 

IV 

1 

120 

1000 

8 

31 

These  tests 
conditions  were 

should  be  considered  like  a pilot  study, 
not  properly  counterbalanced  to  account  for 

The  test 
learning 

and  there  should  have  been  a larger  sample  of  subjects  and  repetitions 
before  making  any  judgments  about  differences  between  test  conditions. 

For  example,  the  RMS  tracking  error  was  recorded  during  the  tests  but 
there  was  no  significant  difference  between  conditions.  In  fact,  they 
were  all  about  the  same. 

On  the  other  hand,  the  plotted  results  on  Figures  9 and  10  seem  to 
validate  the  assumed  human  transfer  function.  The  gain  tends  to  be  con- 
stant with  frequency,  i.e.,  within  the  90  percent  band  of  2 db  except 
for  the  last  couple  of  frequency  bands  where  the  measurement  energy  might 
be  too  low. 

The  phase  tends  to  fall  off  linearly  with  frequency.  Actually,  the 
phase  follows  the  Padfe  approximation,  shown  by  a dashed  line,  more 
closely  than  a pure  delay  which  would  have  slope  of  -0.3  all  the  way. 
Another  observation  is  that  the  assumed  transfer  function  is  a better 
fit  for  the  test  conditions  with  stiffer  legs,  i.e.,  = 500,  1000, 

and  °>. 
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VIII.  DISCUSSION 


Looking  back  over  what  has  been  done,  I conclude  that  there  is  merit 
to  the  analysis  of  this  control  problem.  The  exact  human  transfer  func- 
tion is  probably  very  similar  to  the  one  used  in  these  analyses.  The 
optimum  parameter  values  are  probably  within  the  range  indicated  here, 
i.e.,  K_/B  between  8 and  14  and  /JCp/J  equal  to  30  or  more.  Moreover, 
it  is  possible  to  analyze  such  a problem, then  confirm  the  results  by 
simulation.  Eventually  these  results  should  be  supported  by  field 
experiments  with  hardware. 

The  area  deserves  more  attention  than  it  has  received.  There  are 
any  number  of  applications  in  addition  to  the  military  pointing  task 
which  inspired  the  current  effort.  The  motion  picture  tripod  is  the 
most  similar  task  but  there  are  other  tasks  such  as  remote  manipulation 
of  a mass  through  mechanical  linkage.  This  problem  is  a little  different 
because  the  compliance  is  between  the  point  of  application  of  torque  and 
the  load  rather  than  between  the  load  and  ground.  Where  should  the  vis- 
cous damping  be  introduced  to  stabilize  this  loop? 

There  are  many  things  to  do  to  continue  the  work  started  here. 

Some  of  them  would  improve  the  quality  of  results.  Some  of  them  would 
explore  the  simplifying  assumptions  made  in  the  introduction.  Others 
would  expand  the  generality  of  the  results. 

a.  The  torque  input  device  could  be  built  to  simulate  a specific 
piece  of  hardware  so  that  there  could  be  field  test  data  available  for 
comparison.  The  actual  torques  and  visual  angles  could  be  duplicated 
in  the  simulation.  The  reference  task  would  need  to  be  duplicated  as 
well. 


b.  Once  a good  experimental  setup  has  been  established  there 
should  be  more  data  and  a more  systematic  experimental  procedure  should 
be  followed. 

c.  The  results  should  be  expanded  to  include  other  forcing  func- 
tions. Will  the  transfer  function  and  the  optimum  design  change  much 
if  the  task  is  changed  from  white  noise  to  step  functions  or  to  low 
frequency  sine  waves? 

d.  Other  methods  of  introducing  torque  could  be  considered,  e.g., 
two  hands  as  with  a steering  wheel  and  twisting  as  when  removing  a jar 
lid. 


e.  The  limits  of  linearity  of  the  human  transfer  function  could  be 
defined  both  with  respect  to  visual  threshold  and  maximum  torque. 

f.  Coulomb  friction  should  be  evaluated.  It  is  expensive  to  mini- 
mize coulomb  friction  in  a device  so  we  should  know  just  how  much  can 
be  tolerated  when  there  is  a given  level  of  viscous  friction. 
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A human  operator  is  very  adaptable  and  he  can  usually  manage  to  get 
by  with  a less  than  optimum  tracking  device.  But  a near  optimal  design 
will  require  less  training;  it  will  be  easier  to  use;  and  it  will  give 
better  performance  when  it  really  counts.  Pursuing  the  research  described 
above  will  help  achieve  better  designs  and  perhaps  even  save  money  by 
preventing  over  design  when  it  does  not  help. 
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